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Voltage-gated ion channels are crucial for electrical activity and chemical signaling in a variety of cell types. Struc- 
ture-activity studies involving electrophysiological characterization of mutants are widely used and allow us to 
quickly realize the energetic effects of a mutation by measuring macroscopic currents and fitting the observed volt- 
age dependence of conductance to a Boltzmann equation. However, such an approach is somewhat limiting, prin- 
cipally because of the inherent assumption that the channel activation is a two-state process. In this analysis, we 
show that the area delineated by the gating charge displacement curve and its ordinate axis is related to the free 
energy of activation of a voltage-gated ion channel. We derive a parameter, the median voltage of charge transfer 
(V m ), which is proportional to this area, and prove that the chemical component of free energy change of a system 
can be obtained from the knowledge of V m and the maximum number of charges transferred. Our method is not 
constrained by the number or connectivity of intermediate states and is applicable to instances in which the ob- 
served responses show a multiphasic behavior. We consider various models of ion channel gating with voltage- 
dependent steps, latent charge movement, inactivation, etc. and discuss the applicability of this approach in each 
case. Notably, our method estimates a net free energy change of approximately —14 kcal/ mol associated with the 
full-scale activation of the Shaker potassium channel, in contrast to —2 to — 3 kcal/mol estimated from a single 
Boltzmann fit. Our estimate of the net free energy change in the system is consistent with those derived from de- 
tailed kinetic models (Zagotta et al. 1994./. Gen. Physiol. doi:10.1085/jgp.l03.2.321). The median voltage method 
can reliably quantify the magnitude of free energy change associated with activation of a voltage-dependent system 
from macroscopic equilibrium measurements. This will be particularly useful in scanning mutagenesis experiments. 



INTRODUCTION 

Voltage-gated ion channels are integral membrane pro- 
teins that, upon sensing a change in the membrane 
electric field, open a passage for ions to flux through 
the membrane. They are represented in all major king- 
doms of life and are crucial for both electrical and chemi- 
cal signaling pathways in higher organisms (Hille, 2001). 
Several inherited diseases such as arrhythmias and epi- 
lepsies have been shown to be correlated to mutations 
in these proteins, thereby underscoring their physiologi- 
cal importance (Lehmann-Horn andJurkat-Rott, 1999) . 
To understand the mechanisms of ion channel gating 
and function, it is necessary to obtain accurate estimates 
of energetic effects of site-specific mutations. One widely 
used approach involves measurements of macroscopic 
ionic currents for a series of voltage steps from which 
one can derive the relative fraction of open channel at 
various potentials (Po-V) . These responses show a sig- 
moidal voltage dependence and are typically character- 
ized by fitting to a single Boltzmann equation. For such 
a curve, the chemical free energy difference between 
the open and the closed state is defined by AG C , which is 
characterized by two parameters, z, the Boltzmann slope, 
and V1/2, which is the voltage that elicits the half-maximal 
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response. This free energy difference is also referred to 
as the free energy of channel opening at zero voltage. 

However, the voltage-gated ion channels are known to 
transit through several intermediates to reach their final 
open states (Cole and Moore, 1960; Vandenberg and 
Bezanilla, 1991; Bezanilla et al., 1994; Zagotta etal., 1994a; 
Schoppa and Sigworth, 1998). The slope factor, which 
for a two-state process is the charge translocated, in the 
case of a multistate process depends not only on the 
charge displaced during activation but also on stabilities 
of the intermediate states (McCormack et al., 1991; 
Sigworth, 1994; Bezanilla, 2000; Villalba-Galea etal., 2008). 
Thus, mutations that affect only the energies of the inter- 
mediate states would generate erroneous values of the 
net free energy difference between the initial closed and 
the final open states when derived from Boltzmann fits. 

Furthermore, when the central pore opens, the channel 
molecule has already undergone significant conformational 
transitions as indicated by separation in the charge-voltage 
(QV) and conductance-voltage (GV) curves (Armstrong 
and Bezanilla, 1973; Perozo et al., 1993) . Therefore, the 
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free energy associated with the activation process cannot 
be fully captured from conductance-voltage relationships. 
Single Boltzmann fits of GV curves for Shaker potassium 
ion channel give an estimate of —2 to —3 kcal/mol in the 
chemical free energy difference between the open and 
closed state (Monks et al., 1999; Li-Smerin et al., 2000a; 
DeCaen et al., 2008). This estimate, which corresponds 
to a single hydrogen bond or salt bridge interaction, is 
surprising and seems to suggest that a single interaction 
site could tip the balance between the two states. 

Finally, calculating the free energy change using single 
Boltzmann fits becomes even more problematic when 
mutations result in discernibly multiphasic Po-V curves 
(Li-Smerin et al., 2000b; Yifrach and MacKinnon, 2002). 
Boltzmann fits to the Po-V curve also neglect the energetic 
effect of nonunity values of the maximum open proba- 
bility, which has been observed in voltage-dependent ion 
channels (Schoppa et al., 1992). These observations 
prompted us to carefully analyze the free energy princi- 
ples associated with a voltage-dependent process at a theo- 
retical level. Our objective was to derive a more general, 
physically consistent measure of free energy change that is 
easily obtainable from experimental observables. 

In this paper, we analyzed the fundamental basis of 
free energy change in a voltage-dependent system from 
the standpoint of classical statistical mechanics. The pio- 
neering work of Wyman and Gill (1990) on theoretical 
thermodynamics has greatly furthered our understanding 
of ligand binding equilibria. Of specific interest to us was 
the concept of median ligand activity. The elegant deriva- 
tion of Wyman (1964, 1967) shows that from the knowl- 
edge of the maximum number of ligands that bind a 
macromolecule and the median ligand activity, one can 
obtain an accurate estimate of the mean ligand binding 
affinity. Here, we adopted a similar approach to analyze 
the free energy relationships of voltage-dependent ion 
channels. We found that a comparable parameter, the 
median voltage of activation, as calculated from the QV 
curve, can be related to the net chemical free energy 
change associated with voltage-dependent activation of 
these proteins. Specifically, this procedure yielded a 
value of —14.6 kcal/mol associated with full-scale chan- 
nel activation of the Shaker potassium channels and 
— 16.1 kcal/mol for the voltage-gated sodium channels. 
We have discussed several characteristic models of ion 
channel activation such as those with inactivation, coop- 
erativity, voltage-independent steps, latent charge move- 
ment between multiple open states, etc. and illustrate the 
validity and robustness of our proposed method. 

MATERIALS AND METHODS 

Oocyte expression 

The cDNA voltage-gated Shaker potassium channel, with a dele- 
tion in the N-terminal region (residues 6-46) and with the W434F 
mutation, was linearized with NotI (New England Biolabs, Inc.) 



and in vitro transcribed to generate the mRNA using the mMES- 
SAGE mMACHINE T7 kit (Life Technologies). Stage V Xenopus 
laevis oocytes were injected with 50.6 nl mRNA (at a concentration 
of 0.1 ug/ul). After injection, the oocytes were kept at 18°C in a 
solution containing 100 mM NaCl, 2 mM KC1, 1.8 mM CaCl 2 , 
1 mM MgCl 2 , 5 mM Hepes, 0.1 mM DTT, and 0.2 mM EDTA, sup- 
plemented with 100 ug/ml gentamicin and 100 mg/ml bovine se- 
rum albumin. Measurements were performed 1-2 d after injection. 

cDNAs of both the a and fi subunit of the rNa v 1.4 were tran- 
scribed as described in the previous paragraph. Equimolar ratios 
of the a and (3 subunit mRNAs were coinjected into Xenopus oo- 
cytes to a final volume of 50 nl. Injected oocytes were preserved as 
in the previous paragraph, and measurements were performed 
3-5 d after injection. 

Gating current measurements 

The gating current measurements were performed on a cut-open 
oocyte voltage clamp set-up (CA-1B; Dagan Corporation) as de- 
scribed previously (Muroi et al., 2010; Lacroix and Bezanilla, 
2011). For the potassium channel gating current measurement, 
the external solution was 115 mM NMG-MES (JV-methyl-D-glu- 
camine methanesulfonate), 2 mM Ca-MES, and 10 mM Hepes, pH 
7.4. For the sodium channel gating current measurement, the ex- 
ternal solution was 115 mM Na-MES, 2 mM Ca-MES, and 10 mM 
Hepes, pH 7.4. In the latter case, all ionic currents were blocked 
by the application of 10 pM tetrodotoxin to the external and middle 
chambers. For both channels, the internal solution was 115 mM 
NMG-MES, 2 mM EGTA, and 10 mM Hepes, pH 7.4. The record- 
ing pipette resistance was 0.3-0.5 Mfi. Analogue signals were sam- 
pled at 250 kHz with a Digidata 1440 interface (Molecular Devices) 
and low-pass filtered at 10 kHz. The capacitive transient currents 
were subtracted online using the P/4 method with a subtraction 
holding potential of — 1 20 mV for the potassium channels and 
50 mV for the sodium channels. Gating currents were obtained by 
applying a depolarizing pulse (50 ms for potassium and 20 ms for 
sodium channels) to voltages from — 120 to 10 mV (at 5-mV inter- 
vals) for the potassium channels and —160 to 30 mV (at 10-mV 
intervals) for the sodium channels. The holding potential was 
—90 mV, and a 50-ms-long pre- and postpulse at —130 mV was used. 

Analysis and simulations 

Each of the QV curves is a mean obtained from five oocytes. The 
curve was fitted to various equations as described in the Table 1. 
Computation of the area was performed using the trapezoid 
method. In brief, for a fractional gating charge displacement ver- 
sus voltage (QV) curve, with n points, the area between the curve 
and the ordinate (Q) axis was calculated as 

where V is the f" point on the QV curve and Q is the fraction of 
charge that is transferred at voltage V,. Numerical simulations of 
the different models were performed using MATLAB 2008b 
(MathWorks) . 

RESULTS 

Theory 

Free energy change in a two-state process is obtained from 
a Boltzmann ■fit. Consider a voltage-dependent ion chan- 
nel that can exist in two states, closed (C) and open (O). 
On increasing the voltage, the open state of the channel 
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becomes increasingly populated. The associated free 
energy difference between the closed and open states can 
be written as AG C - qV, where V is the membrane elec- 
tric field gradient and AG C is the chemical (nonelectri- 
cal) free energy difference between the two states or the 
free energy change associated with the transition in the 
absence of an electrical driving force (V = 0) . q is the gat- 
ing charge translocated when the ion channel activates 
and is responsible for an electrical component in the 
net free energy change associated with activation. The 
equilibrium constant for such a transition at any voltage 
will be 

^ = X°exp( 9 Vp), 
"c 

where P 0 and P c are the fractional occupancy of the 
open and closed states, respectively, and K° is the voltage- 
independent component of the equilibrium constant, 
which is related to AG C as AG C = -RT\nK°. The voltage- 
dependent probability of occupancy of the open state, 
in this situation, will be 

P - P o - 1 1 

° P 0 + P c 1 + {P C /P 0 ) exp(- ? Vp) - 

K° 



Rewriting K° as exp( — gfVi/g), where V1/2 (Stevens, 1978) 
is the voltage at which half of the channels are open, P 0 
can be reexpressed as 



0 1+ exp{q(V 1/2 - V)pf 

By fitting an experimentally derived Po-V curve to Eq. 1 , 
one can derive the parameters of the Boltzmann fit, 
V1/9, and (/and thereby estimate AG C as qFVi/^. A similar 
equation can be derived for the gating charge displace- 
ment versus V relation. If a channel can access only two 
states, the normalized QV and Po-V curves will super- 
impose on each other. Also, the Boltzmann slope of the 
curves will equal the exact number of gating charges 
displaced during channel activation. However, voltage- 
gated ion channels are known to transit through a large 
number of intermediates during activation; the experi- 
mentally derived Po-V and QV curves do not superpose, 
with the latter almost always preceding the Po-V curve on 
the voltage axis. Also the Boltzmann slope of either curve 
is much lower than the actual amount of gating charge 
displaced during channel activation (Schoppa et al., 1992; 
Aggarwal and MacKinnon, 1996; Seoh et al., 1996) . Thus, it 
is likely that free energy estimates of ion channel activation 
via a Boltzmann fit to activation response are inaccurate. 



TABLE 1 

Estimation of free energy change from QV 
Fitting function Parameters Sh-IRW434F Na v 1.4 

Values . ~ Values . ~ 



1 

l + exp)zF(V 1/a -y)PI 



[l + exp{zF(V 1/2 - V)P)] 4 



(z,+ z 2 )^l + explz^V™- V)P) 1 + exptz^fV*'- 



z 1 exp|z 1 F(y i ' / 1 »-l/)PI + ( h +z 2 )exp[{z 1 V 1 ' / 1 j + h V™ - ( h + z, )V }F$] 



Zj + z 2 



1+ exp{z 1J F(y i ( / 1 »-y)P} + exp^y™ 



? M/2 



(z, + z 2 )y}^p] 



F]^* VdQ = Q ma FV Mt 



Vi/s 


-45.9 


-3.2 


-56.5 








z 


3.03 


1.4 


v 1/2 


-64.6 


-13.93 


-97.15 




2.35 




1.08 


y 1/2 (1) 


-58.2 




-71.6 


Zl 


1.99 


-7.07 


1.1 








-42.4 


-49.64 




4.5 




1.85 


y 1/2 (1) 


-56.64 




-71.02 


Zl 


2.02 


-6.44 


1.15 


y ]/2 (2) 




-42.6 


-48.23 


-2 


3.89 




1.58 


v„ 


-46.81 


-14.63 


-56.28 










13.6 


12.8 



-1.82 



-9.63 



-3.92 



-3.62 



-16.5 



The table lists the different functions used to fit the QV curve and the values of the corresponding parameters derived from the fits to the QV curves of 
the Sh-IRW434F and Na v 1.4 channels. V1/2 parameters have units in millivolts, and 2 parameters have units of electronic charges. The free energy changes 
calculated from each of the fits are also listed (kcal/mol). 
a A two-state model and the fitting function is a Boltzmann function. 

b A model where the channel activates in four independent single-step transitions; the fitting function is the fourth power of the Boltzmann function and 
AG C =4zFV V2 . 

C A model where the channel activates via two energetically independent steps; the fitting function is the sum of two Boltzmann terms and 
AG C =z I fy™ + h FV$. 

d A three-state model where the channel activates in two sequential steps; AG C = z^FV^l + z^FVy^, 

e The free energy change calculated via the median method of channel activation. The Qr, ax values for each channel was obtained from previously published 
literature (Schoppa et al., 1992; Hirschberg et al., 1995; Aggarwal and MacKinnon, 1996; Seoh et al., 1996). 
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Figure 1 . Free energy measures from a QV curve. (A) A hypo- 
thetical normalized QV curve. The vertical dashedjine is^the V = 0 
axis. The shaded area is the integral J VdQ, ( Qy = Q / <2m ax ), 
which equals AG C / Q^^. The integral is negative on the left of 
the V = 0 axis and positive on the right, as indicated by the — and 
+ signs. (B) The same QV curve with the median voltage of charge 
movement (V m ) marked by the vertical dashed lines. The areas on 
either side of the V m axis are equal and are both positive. (C) The 
QV curve intersected by the V = 0 and V = V m axes. The area of the 



Area under the activation curve is a measure of the free en- 
ergy change of activation. Consider that the voltage-gated 
ion channel exists in n different conformational states, 
with each state i associated with a gating charge (or va- 
lence), given by 



(2) 



where the summation is over all the charges z, of the 
channel protein in state i, which are informed princi- 
pally by the charges on the side chains of specific ioniz- 
able residues of the protein. X,, is a state-dependent 
parameter representing the fraction of the membrane 
electric field sensed by the f h charge of the protein in 
the state i (Stevens, 1978; Roux, 1997). q b as described 
by Eq. 2, essentially governs the voltage dependence of 
the energy of a given conformational state, E„ i.e., 



3^ 
dV 



(3) 



The negative sign in Eq. 3 implies that states with a 
more positive valence become more stable with increas- 
ing voltage. 

The mean gating charge of the entire ensemble of con- 
formational states at any given voltage is (Sigworth, 1994) 



(4) 



where P,(V) is the equilibrium occupancy of the i state 
at voltage V, which is assumed to follow a Boltzmann 
distribution. Thus, P t = exp(-E$)/Z, where (3 = l/&sT 
and Z is the partition function of the system, 



Z = ^exp(-Efi). 



The mean free energy of the ensemble, assuming that 
temperature is held constant, is 



k B T\nZ. 



(5) 



Changing the electrical potential (V) will alter the 
distribution of the channel states and thus the mean 
ensemble free energy. At a constant temperature, dif- 
ferentiating Eq. 5 with respect to voltage gives 



3G 
dV 



(6) 



shaded rectangle is equal to the sum of the two areas shaded in A, 
taking their respective signs into consideration. 
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Now using Eqs. 3 and 4, the differential change in the 
free energy of the system (Eq. 6) with voltage can be 
expressed as 



dv * 



(7) 



Thus, network done or net change in free energy of the 
system as the voltage is changed from Vi to V2 can be 
given by 



AG= G 



rV 2 _ 

- QdV. 



(8) 



The gating charge displacement measured when the 
voltage is switched from a reference voltage (V re f) to a 
voltage V is essentially the difference Q(V) - Q(V ref ). V re f 
is usually a hyperpolarizing voltage when all protein 
charges are retracted to their initial resting configura- 
tion, and, without any loss of generality, Q(V„ f ) can be 
taken to be 0. Taking Vi as V re f and using V instead of V?, 
Eq. 8 can be integrated by parts, giving 



AG-- 



-(QV) + 

' v ref 



•Q(V) 



VdQ. 

J Q.(v ref ) 



('.)) 



Suppose gating charge movement is measured experi- 
mentally from a hyperpolarizing reference potential, 
where all the charges are retracted, to a potential V at 
which all charge movement saturates. Thus Q(V) at 
highly depolarized voltages is Qma*, which makes Eq. 9 



AG = -Q max V + 



* Qma: 



VdQ. 



(10) 



From Eq. 10, the net free energy change in the process 
thus comprises two terms: the first (— QmaxV) scales lin- 
early with voltage, and the second component ( AG C ) 
equals the area between the gating charge displace- 
ment versus voltage (QV) curve and the ordinate (Q) 
axis. A hypothetical QV curve for a voltage-dependent 
ion channel is shown in Fig. 1 A, where the shaded re- 
gion represents AG C . The linear component is the 
change in the net electrical energy of the system as it 
undergoes a complete transformation. The second com- 
ponent, unlike its linear counterpart, is convergent and 
will be shown to lead to an expression for the net change 
in the chemical free energy of the system. In the next 
section, we derive a parameter, the median voltage of 
activation, which is directly related to AG C . 

The median voltage of activation. According to the defini- 
tion of median ligand activity proposed by Wyman 
(1967), we define the median voltage of activation, V m , 
as the voltage at which 



(On 



Q)dV. 



Also let Qm be the gating charge displacement at V m . 
By definition, AG C is 

AG C = \ VdQ = \ (V - VJdQ + Q mzyi V m . (12) 
Jo Jo 



The integral on the right side of Eq. 12 can be broken 
down into two integrals: 



(13) 



ra™, - r&n - 

I (V-VJdQ = I (V-VJdQ 
Jo Jo 

- P m " (y-vjdfo^-o). 

J £>... 



Integrating each of the two integrals in Eq. 13, by parts, 
and using the definition of V m (Eq. 11) and the fact that 
Q(K») = Qm> it can be shown that the integral on the 
right of Eq. 12 is 0. The result is that 



(14) 



Fig. 1 B shows the QV curve intersected by the V = V m 
axis, with the two shaded regions being equal in area. 
The areas on the right and left of the V m axis correspond 
to the integrals on the right and left side of Eq. 11, respec- 
tively. The area of the rectangle shaded in Fig. 1 C is equal 
to the sum of the two areas shaded in Fig. 1 A, taking the 
sign of the latter areas into consideration. 

A physical interpretation of AG C and V m . Thus far, we have 
shown that the area under an experimentally measured 
QV curve can be separated into two components, a lin- 
ear and a saturating (converging) component, both be- 
ing dependent on the total charge moved between the 
initial and final states of the system. The saturating com- 
ponent, AG C , can be described (Eq. 14) by a parameter, 
the median voltage of activation, given by Eq. 11. In this 
section we will show that AG C is the change in chemical 
free energy of the system associated with the activation 
of the system. 

According to Eq. 3, the energy of each state of the 
system can be written as E t = E? - q { V, where E? is the 
chemical energy of the 2 th state. The Boltzmann weight 
of each of these states can be expressed as K { exp(^ VP) , 
where Kj = exp ( —E$) . We assume that at sufficiendy hy- 
perpolarizing voltages, the system occupies just one 
(reference) state for which q r = 0 and E r = 0. With this 
normalization, the Boltzmann weights of the different 
states become exp (?,VP), where K° r ^ (exp(-|i , I -£,°||3)) 
is the chemical component of the equilibrium constant 
of the transition from the reference state to state i. The 
partition function of the system is now 



(11) 



Z = exp ( ?i Vp). 



(15) 
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For our reference state, K° T = 1. 

Now we turn to Eq. 1 1 and rewrite it as 



n o.dv = r q„ 



,dV. 



Combining Eqs. 5, 7, and 16, we get 
3lnZ 



k > T L ~^r dV 



or 



k„TlnZ\ 



(16) 



(17) 



(18) 



At hyperpolarizing voltages (V — > — °°), only the refer- 
ence state will be populated, and thus lim^^lnZ = 0. 
Dividing the equation by k s T and rewriting QmaxV(3 as In 
ex p(QmaxV(3), where (fegT = 1/(3), Eq. 18 becomes 



lim T , 



In 



Z 



V 

max ' m 



(19) 



Now, we substitute Z from Eq. 15 into Eq. 19. At very 
high voltages, all the terms of Z are dominated by the 
(final) state. Thus, the integral equations, Eqs. 16 and 
17, ultimately yield 



Q V 



■k R T\nK° 



(20) 



where K r ^ f is the chemical component of the equilib- 
rium constant for the transition from the reference 
state to the final state of the system. Two fundamentally 
crucial points emerge from the relation described in 
Eq. 20. First, AG C , which is equal to QmaxV m , is thus the 
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Figure 2. QV curves of the Shaker potassium channel and the 
rNa v 1.4 channel. Normalized QV curves for Sh-IR W434F and 
Na v 1.4. Each curve is a mean of independent measurements from 
five oocytes. The normalized QV curve for the ILT mutant was 
simulated according to the kinetic model proposed by Ledwell 
andAldrich (1999). 



measure of chemical free energy difference between the 
initial resting state of the system and the final activated 
state of the system, i.e., the free energy difference be- 
tween the initial and final states of the system in the 
absence of any electric field (0 mV) . Second, at V m , the 
Boltzmann weight of the fully activated state of the system, 
X r °_^exp(Q max V„P), is 1, which is same as the Boltzmann 
weight of the reference state (for which the valence is 0) . 
This implies that at V m , the state or states where all 
charge has moved (fully activated) and that where none 
has moved (fully resting state) are equally populated. 

The relations established are independent of the 
pathway connecting the initial and final states of the sys- 
tem and the number of possible states of the system. 
Measurement of the QV can thus directiy be related to 
the change in the free energy without elaborate model 
fitting procedures or assumptions about the nature 
of the conformational transitions (such as the two- 
state process assumption required by Boltzmann fits) . 
The relations remain valid even when the QV curves 
show biphasic or asymmetric behavior. It requires the 
knowledge of Q mx , which, at least for the voltage- 
gated potassium (K v ) channels, is known with reason- 
able certainty (Schoppa et al., 1992; Aggarwal and 
MacKinnon, 1996; Seoh et al, 1996). Thus, at a first 
level of approximation, the measurement of the me- 
dian voltage of activation would facilitate a straight- 
forward and accurate calculation of the free energy 
changes associated with the voltage-dependent activa- 
tion of any system. 

Accuracy of Q max V m as a measure of free energy change 

To compare the free energy values derived from QmaxVm 
and the Boltzmann fits for real systems, we experimen- 
tally recorded the gating currents in the inactivation- 
removed Shaker W434F potassium channel construct 
(Sh-IR W434F) and in the voltage-gated sodium chan- 
nel (Nayl.4), with the ionic currents blocked with tetro- 
dotoxin, and evaluated their normalized QV curves 
(Fig. 2). V1/2 for each was obtained by fitting the 
Boltzmann equation (Eq. 1), whereas V m was obtained 
by calculating the area between the curve and the 
Q axis (see Materials and methods). Previously pub- 
lished values of Qmax for each channel (Schoppa et al., 
1992; Hirschberg et al., 1995; Aggarwal and MacKinnon, 
1996; Seoh et al., 1996) were used for the free energy 
calculations using the median method. We find that the 
free energy values computed via the median method 
were four to eight times the values computed via the 
Boltzmann method (Table 1). 

Occasionally for potassium channels, the sigmoid activa- 
tion curves have been fitted to a fourth-power Boltzmann 
equation (Zagotta et al., 1994b; Smith-Maxwell et al., 
1998b; Gonzalez et al., 2000, 2001), assuming that iden- 
tical conformational changes occur in a single step in 
each of the four identical voltage-sensors, which "move" 
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TABLE 2 

Comparison of free energy estimates from kinetic models and median method 



Channel 


X? 


Zl 


K 


28 


K 




AGc = - RT In { ( K° K° ) 4 K° L | 




v„, 


Q„ la xFV m 
















kcal/mot 




mV 


kcal/mol 


Sh-IR 


3 


1.8 


132.07 


1.4 


12 


0.4 


-15.58 


13.2 


-46.8 


-14.19 


ILT 


74.667 


1.7 


466.667 


1.3 


0.0143 


1.8 


-22.17 


13.8 


-67.4 


-21.37 



The table lists the parameters of the ZHA model for the Sh-IR channel and the ILT mutant. V m value for the ILT mutant was calculated from the QV curve, 
generated using the model parameters published previously (Ledwell and Aldrich, 1999; Fig. 2). The values used in both cases are those specified 
in the ZHA model. 



independent of the other. This fit provides a much 
closer value of the free energy change to that obtained 
from the median method for Sh-IR W434F, but for 
Na v l-4 the difference between the two estimates con- 
tinues to be large (Table 1). Fitting the QV curve of 
Na v l -4 to the product of four nonidentical Boltzmann 
terms (which takes into consideration the four non- 
identical domains of the protein) would be ill constrained. 
We fitted two other model dependent sigmoidal func- 
tions (Perozo et al., 1994; Rodriguez et al., 1998) , both of 
which resulted in free energy estimates significantly 
different from the median method for Sh-IR W434F 
as well as Na v 1.4 (Table 1) . 

For the Shaker potassium channel, elaborate multi- 
step kinetic models have been proposed by various 
groups (Bezanilla et al., 1994; Zagotta et al., 1994a; 
Schoppa and Sigworth, 1998) based on single channel 
conductances, gating currents, and macroscopic ionic 



currents. Here, we consider the model proposed by 
Zagotta etal. (1994a; ZHA model). 



- R 2 ^=^ A 



R 1 s R 2 s A 

K° ; z, K° ; z 2 

R.p R 2 t A 

K° ; z, K° : z 2 

R-p R 2 ^ A 



C O 



(SCHEME 1) 



According to this model, each voltage sensor inde- 
pendently activates in two discrete steps (with equilibrium 
constants K\ and K?) , and when all the voltage-sensors 




Figure 3. Effect of an intermediate voltage-independent step on the free energy estimates. (A) A six-state linear model of activation of a 
voltage-dependent ion channel with five closed state and a single final open state. The voltage dependence of the equilibrium constants 
is given by K i = Kf exp(z, FV$) (i = 1, 2, 3, 4, or 5), where K° is the chemical component of the equilibrium constant K, . and z, is its 
voltage dependence. The third step (shown in bold) is voltage independent, and thus z 3 = 0. AG C is the sum of the chemical free energy 
change of each of the steps and will equal — .RTlnJ^J X?. (B) QmaxFV m is plotted against AG C for different values of log (—5 to 5) 
and two values of z% (closed squares, z% = 5; and open triangles, z% = 0) . The remaining parameters of the model were arbitrarily chosen as 
K° =1, K° = 75, Kl = 25, K° = 5, z x = 2, z 3 = 0, z 4 = 2.5, and z 5 = 1. (C and D) Variation of AG C (C) and V m (D) with changing values 
ofz;, (0-5) and log K°. 
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are activated, a concerted transition gates the ion channel 
pore open (equilibrium constant K L ). All the equilib- 
rium constants have an exponential voltage dependence: 
K { = Kfexpi^FV^) (i = 1, 2, or L). The net chemical 
free energy change in going from the reference state to 
the final open state, AG C , will be -RTln{(K° K°f K° L }. 
(Note the use of R, the universal gas constant, instead of 
k B and the inclusion of Faraday constant [F] in the nu- 
merical free energy calculations.) Using the parameters 
reported, we find that the free energy estimate from the 
ZHA model parameters and the median method are 
nearly identical (Table 2) . Next, we considered the ILT 
mutant of the Shaker potassium channel (Smith-Max- 
well et al., 1998a; Ledwell and Aldrich, 1999), which is 
characteristic in that it results in a QV curve with a sec- 
ondary phase, widely separated from the primary (Fig. 2) , 
and thus cannot be fitted to a Boltzmann curve. How- 
ever, the value of chemical free energy change for this 
mutant computed using the median method from the 
numerically simulated QV curve is again almost the 
same as that obtained from the kinetic model parame- 
ters. The fitting exercise described herein thus shows that 
the median method is a physically consistent, experi- 



mentally feasible method to accurately estimate voltage- 
dependent energetics in the system. 

Free energy estimation using the median method in the 
case where certain steps are voltage independent 

Next we considered models of ion channel activation 
with voltage-independent steps to illustrate how our in- 
ferred energetics can be modified by consideration of 
such transitions. Let us assume that a voltage-independent 
transition occurs in the middle of a linear transition 
pathway, involving six states with the final state being 
the open state O t (Fig. 3 A, scheme I) . Because all of the 
steps beyond it are voltage driven, as we increase the 
voltage, the system will be pushed further and further to 
the right until it finally saturates in Of. For this system, 
the chemical free energy difference between the initial 
and final states, AG C , will be equal to QmaxVn as de- 
scribed previously (see Theory section A physical inter- 
pretation of AG C and V m ) . 

To test this numerically, we generated a large number 
of QV curves, based on this model, for different values 
of Z2 and K°. The former is the voltage dependence of 
the equilibrium constant for the second step, and the 
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Figure 4. Effect of terminal voltage-independent steps on the free energy estimates. (A) Two six-state linear schemes of voltage- 
dependent activation of a channel. In both schemes, except for the final transition (in bold) , all other steps are voltage dependent. 
Scheme II has one open state, whereas scheme II has two open states. The voltage dependence of each of the equilibrium parameters 
in schemes II and III is defined similar to those in scheme I (Fig. 3 A) . For schemes II and III, z 5 = 0. (B) Results of the numerical 
simulations performed using scheme II. = 10, Zj = 1.5, 23 = 1, and z 5 = 0. The remaining parameters were same as those in scheme I 
(Fig. 3 A). K 5 was varied from 0.001 to 2,097.15. The relation between the free energy difference between the initial and final states, 



AG C ( -RTlri[i i=1 K? ), and the free energy change in the ensemble, QniaxFV m (closed symbols) and 
different values of K®. Adding Q /ma[ FV m and the correction factor, -RT In 



(open symbols), for 



gives AG C , as depicted by the dashed line. (C) P 0 



values calculated at different values of K 5 for scheme II (A) . 
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latter is the equilibrium constant for voltage-indepen- 
dent third step. For this scheme CL™ = X, z ;> and thus 
changing z 2 alters the maximum gating charge trans- 
ferred in the process of activation. From each simulated 
curve, we computed the median voltage of activation, V m , 
by integrating the area between the curve and the ordinate 
(Q) axis (jl VdQ f ). Plotting Q^FV™ against AG C (Fig. 3 B) 
shows that they are identical to each other and indepen- 
dent of Qmax. The surface plot of AG C for varying z 2 and 
log k" (Fig. 3 C) shows that it has logarithmic depen- 
dence on K° and is independent of z%. However, V m 



shows a logarithmic dependence on iq but an inverse 
dependence on Zg (Fig. 3 D). AG C , which is the chemi- 
cal free energy difference between the final and initial 
states of the system, will be independent of the magnitude 
of gating charge translocated between the two states. 
For a fixed value of Qmax in the system, V m is linearly re- 
lated to AG C , and both will have a proportional depen- 
dence on the equilibrium constants. When Q^a* is altered 
without changing any of the equilibrium constants, AG C 
remains unchanged, but V m will be altered because 
Q^a X FV m = AG C . " 




1 1 ' 1 ■ 1 •— H 1 1 1 ■ 1 1 1 1 1 1 1 

-40 -30 -20 -70 -60 -50 -40 -30 -20 

AG C (kJ/mole) AG C (kJ/mole) 

Figure 5. Influence of latent charge movement and cooperativity on inferred energetics. (A) A 10-state MWC cooperative scheme of 
voltage-dependent ion channel activation. The channel comprises four identical voltage-sensing modules and one pore domain, each 
capable of existing in two conformations. States Q and O, differ in the conformational status of the pore domain, whereas the different 
C,s and 0,s differ among each other in the number of activated voltage-sensing modules ( i = 0, 1 , . . . , 4) . Activation of the each of the volt- 
age sensors facilitates the opening transition of the pore and vice versa. The number alongside each state indicates its multiplicity. The 
voltage dependence of the equilibrium constant for activation of the pore and voltage sensors follows the relation K t = Kf exp ( z, F V$) 
(i = V, P) . (B) Plot of QjnaxFVm against AG C for different values of 0 (varied between 5 and 80) and two values of z v (open triangles, z v = 
1; closed squares, z v = 3). The other model parameters were chosen to be Ky = 20, Kp = 10 -5 , and z P = 1.5. Different values of Zy lead 
to different latent charge movement. (C) Plot of Q^axFVm and —RT\nP^ lax _ against AG C for different values of 0 (varied between 5 
and 80) when Zp = 0. The arrow shows the value of 0 and i 3 ^™ beyond which AG C and QmaxFVm deviate. For these simulations, Ky = 20, Kp = 
10~°, and z v = 3. /g™ in each case was M). Adding QmaxFVm and the correction factor, -RT\nP 0 , gj ves AG C , as depicted by the dashed 
line. (D) Plotof Q^axFV,,, and -RT\n{\ - P 0 min ) against A_G C for different values of K° (varied between 10~ 5 and 50) when z P = 0. The 
arrow shows the value of Kp and P™ 11 below which AG C and Q [!lax FV m deviate. For these simulations, K v = 20, 0 = 20, and z v = 3. 
P™"* in each case was ~1. Adding QmaxFV,,, and the correction factor, -i?Tln(l - P™ n ) , gives AG C , as depicted by the dashed line. 
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Next, we modified scheme I (Fig. 3 A) by making the 
last (rather than the middle) step of the pathway voltage 
independent (Fig. 4 A, scheme II) . Now at saturating 
depolarization, both the final states, C5 and Of, will be 
occupied, their relative occupancies being determined 
by the voltage-independent equilibrium constant of 
the last transition ( K b ). Invoking Eq. 19 for scheme II 
(Fig. 4 A) , we see that 



V = 



-h B T lim v 



In 



.exp(e max l/(3) 
,). 



(21) 



Eq. 21 tells us that if the system does not saturate to a 
single state, QmaxVm reports the change in chemical free 
energy in taking the system from its initial reference state 
to a saturating condition, which is an ensemble rather 
than a single state. This brings us to an important point: 
there are two nonequivalent definitions of free energy 
change in the system. One is the change in the free en- 
ergy of the ensemble, whereas the other is the free energy 
difference between the unique final and initial states of 
the system (AG C , which in this case is _ V^n^q^cy ). 
The former is estimated simply via QmaxVm, whereas to 
estimate the latter, we rewrite Eq. 21: 



(22) 



If Kl is not large, the second logarithmic term in Eq. 22 
would have a nonvanishing contribution to the overall 
equation. In this situation, however, will be equal 
to </(!+ Kl). This gives 



AG r 



) V 

will ax m 



k B T\nP™ 



(23) 



Thus, when the final step is voltage independent, the 
knowledge of the maximum open probability, along 
with Qmax an d Vnu would let us estimate the net free 
energy change associated with transferring the channel 
from its initial state Ci to its final state Of. 

We illustrate this point through numerical simula- 
tions. We calculated the V m for several QV curves, based 
on scheme II (Fig. 4 A) , using different values of K° 
and plotted QmaxFVm against AG C (Fig. 4 B). In each 
case, we also calculated P™ x (Fig. 4 C) . Fig. 4 B shows 
that for large values of K°, Q, nax FV m and AG C are identi- 
cal and -RTItlPq** is ^0 (because P^T* is ^1). This is 
expected because under these conditions, the final sat- 
urated state of the ensemble is populated by a single 
state so that the change in the chemical free energy of 
the ensemble and the chemical free energy difference 
between the final and the initial state are equivalent. 
For smaller values of K h , where P™** < 1, these two free 



energy measures are different, and the value of -RT InP™* 
increases with increasing difference between the two 
free energy quantities. 

Next, we consider scheme III, which is a six-state lin- 
ear model with two final open states, the transition be- 
tween them being voltage independent (Fig. 4 A) . For 
this case, the QmaxVm measure is equal to 



0 V 



(24) 



This scheme yields P™** = 1, and, as evident from Eq. 24, 
QmaxVm will be the free energy change of the ensemble 
but not the chemical free energy difference between 
the first and last states (-k B T\aK^ a ). Unlike our pre- 
vious case (terminal voltage-independent step between 
a closed and open state) , here, the two free energy mea- 
sures cannot be distinguished using P™** values. 

Effects of cooperativity and ion channel inactivation 

In the models considered thus far, the median measure 
of the chemical free energy change of the system corre- 
lates well with the model-based estimates. To further 
check and validate these correlations, we performed 
numerical simulations on an allosteric model of ion 
channel activation. Such models have been used with 
success to understand the voltage-dependent gating of 
the BK channels (Horrigan and Aldrich, 1999, 2002) . 
We used an MWC model (Monod et al., 1965) for our 
simulations here. 

Consider scheme TV, a 10-state MWC model of ion 
channel activation as shown in Fig. 5 A. The channel is 
assumed to comprise four identical voltage-sensing 
modules and a single pore domain. K v is the equilib- 
rium constant of activation of each of the four voltage- 
sensing modules, whereas is the equilibrium constant 
of the pore-opening transition. Both of the equilibrium 
constants are assumed to have an exponential voltage 
dependence: K t = K t exp ( z l FV$) , (i = V, P) where z, re- 
flects the voltage dependence of each transition and 
Kj is the chemical component of the equilibrium con- 
stant of each transition. This model is different from 
those discussed previously, in part, because of the pres- 
ence of multiple open states with different valences, 
which will manifest itself as a latent charge movement 
between the open states (Sigg and Bezanilla, 1997). 0 is 
the cooperativity parameter that represents the facilita- 
tion experienced by the pore-opening process caused by 
activation of a voltage sensor and vice versa. 9 is assumed 
to be voltage independent. In scheme TV (Fig. 5 A) , if all 
the steps are voltage dependent, the system will start 
from state Q at hyperpolarizing voltages and end in 
state Of at highly depolarizing voltages. The chemical 
free energy difference between the states Ci and Of, 
AG C , will be given by -RT ln\(Q K° ) 4 K° p ). 

We generated several QV curves for scheme IV (Fig. 5 A) 
with different values of Zy and 0 when z P is not 0. The free 
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energy change calculated from median voltage of acti- 
vation in each case was compared with AG C . Qjnax for 
scheme IV (Fig. 5 A) is 4z v + z P . The plot of QmaxFVm 
against AG C (Fig. 5 B) shows that, as long as the pore- 
opening transition is voltage dependent (z P * 0) , the two 
free energy measures are equal, dependent on 0 but in- 
dependent of Qmax- This illustrates that the presence of 
multiple open states or latent charge movement between 
the open states does not limit the applicability of the me- 
dian method to estimate the chemical free energy change. 

When the pore-opening process is voltage indepen- 
dent, at saturating voltages, both states Cf and Of will be 
populated, their relative proportions being determined 
by the magnitude of the equilibrium constant, K P Q 4 . In 
addition, at hyperpolarizing voltages, both the states Co 
and Oo will be populated, their relative proportions being 



determined by Kp. Thus, at both hyperpolarized and 
depolarized conditions, the system does not exist in a 
unique state but rather in an ensemble of states. From Eq. 
18 and assuming Co to be our reference state, we see that 



<2max^m = " + ( lin V->~ Qmax V) = 

- k B T In ( K° c ^ Cj + K° Cq _ 0/ ) + k B T In (1 + ^ 



(25) 



In terms of the model parameters for scheme IV (Fig. 5 A) , 
Eq. 25 may be rewritten as 



Q***V m = - k B T In { K° P (K° V Q) 4 } - k B T In 



^ + TS° A 4 



(26) 



+k B Tln(l+K° p ). 
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Figure 6. Influence of inactivation on the median estimate of free energy change. (A) An eight-state model for an inactivating chan- 
nel. The encircled states represent the initial reference state (V r P c Ir, red) and the final state (V a PoIa. black) of the system. In the final 
state, all of the units of the system are in their activated/open conformations, whereas in the reference state, all are in their resting/ 
closed conformations. The equilibrium constants are assumed to have an exponential voltage dependence, K i = Kf exp ( z i F Vf5) , (» = 
1, 2,..., 12). (B) An equivalent nested coupled model of inactivation. The conformational transitions in each of the structural units are 
K t = A'° exp ( Zj F V$) , where i is V, P, or inact. The pairwise coupling parameters, Oyp, 8pi, and 9iv, denote the voltage-independent cou- 
pling between units V and P, P and I, and I and V, respectively. The parameters of the cubic model can be related to those of the nested 
model. For instance, K w = i^ nact and K n = Ai„ ac t6pi0iv- (C) Plot of QmaxFVm against AG C for different values of K inacl (varied between 
0.0001 and 0.0256) and two values of Zj nac t. The remaining parameters were arbitrarily chosen to be K v =75, Kp =1, Gyp = 5, 0 P i = 8, 
6iv= 100, z v = 3, and z P = 1. Here, AG C = - RT In K^ Pcls ^ p 1 =- RT In (£° K° K^) = -RTln(Ky Kp K^Qyp Q PI 6^). (D) Maximum 
fraction of inactivated channels at depolarizing voltages for different values of K inacl . 
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For such an MWC model, P™"* at highly depolarized 
voltages will be Kp 9 4 / ( 1 + K° P 0 4 ) , whereas at very low hy- 
perpolarizing voltages, -P™ n will be Kp /(1 + Kp). Using 
these definitions in Eq. 26, we obtain 

AG C = Q Jsm V m - k B T\nP™ + k B T\n(l - 



(27) 



For illustration purposes, we compared numerically 
simulated values of QmaxFV m and AG C in the case when 
the pore opening is voltage independent (i.e., z P = 0) at 
different values of the model parameters. When K P is 
small (thus, io" ln is ^0) , and 0 is large, PS 13 * is ^1 and 
the two free energy measures are equal (Fig. 5 C). For 
smaller values of 0, AG C can be obtained by adding the 
correction factor -RThiPg** to QmaxFVm. Alternatively, 
when both 0 and K P are large (thus, -P™ is significant) , 
the deviation between the two free energy measures can 
be accounted for by the correction factor, i?Tln(l - P™ n ) 
(Fig. 5 D) . Thus, using the median measure along with 
the appropriate correction factors, as described in Eq. 27, 
the chemical free energy difference between the termi- 
nal states can be estimated for systems undergoing volt- 
age-dependent activation according to an MWC scheme 
of activation. 

We finally consider the case of channels that undergo 
inactivation on sustained depolarization (Fig. 6 A, 
scheme V) . The system can be modeled by considering 
that there are three specialized units, V, P, and I, each 
capable of existing in two conformations. Conforma- 
tional change of V (Vr to Va) governs the initial voltage- 
sensing steps of the channel, that of P (Pc to Po) is the 
gate opening transition, and that of I (Ir to Ia) controls 
the inactivation process. Scheme rV (Fig. 6 A) has 12 
parameters. This cubic model of an inactivating chan- 
nel can be represented as a nested coupled model as 



depicted in Fig. 6 B. Each of three structural units has 
an intrinsic activation constant described by a voltage- 
independent component, K i , and voltage-dependent z ; 
(where i is V, P, or inact) . Additionally there are three 
coupling parameters, 0yp, 6pi, and 0rv, which are voltage 
independent and depict the pairwise interactions be- 
tween the structural units. Each of the 12 parameters of 
the cubic model (Fig. 6 A) can be uniquely related to 
the parameters of Fig. 6 B using the principle of micro- 
scopic reversibility. 

If all of the equilibrium constants in this model have 
a positive voltage dependence, at highly depolarizing 
voltages, the system will saturate in the state VaPoIa, 
whereas at very low voltages, only the initial reference 
state, VrPcIr, will be populated. QnaxVm (V m being ob- 
tained from the QV curve) will thus be the chemical 
free energy difference between the initial and final states 
of the system. 

Next, assume that the inactivation process is voltage 
independent. As a result, all of the equilibrium con- 
stants between two connected states differing in the 
conformation of I (Fig. 6 A) will be voltage independent, 
i.e., Kg, Kio, K n , and do not change with voltage. 
Because all other transitions are voltage dependent, both 
states VaPoIr and VaPoIa will be populated at depolar- 
ized potentials, and at very low voltages, both VrPqIr and 
VrPcIawUI exist. This situation is similar to our description 
of the MWC model, with voltage-independent pore open- 
ing transitions, in that multiple states are populated at 
both hyperpolarized and depolarized voltages. 

Thus, Eq. 18 here transforms to 



-k B Tln 



jf° jf° _i_ jf° jf° jf° 



0 \ 



1 + K" 



(I 



V . 

max m 



(28) 



TABLE 3 





Relation between AG C 


and <K ../ \ ,. for ion channel gating models 


Model 


Reference state 
(r) 


Final state 
(f) 


Relation 


Noninactivating; all steps voltage 
dependent or intermediate 


Ci 


o, 


AG C = -flTlnX^ 0/ = Q^FV, 


step voltage independent 








Noninactivating; last step voltage 
independent 


Q 


o, 


AG C = -RTlnK°^ 0f = Q m .^FV m - i?TlnP™" 


Noninactivating; MWC scheme; 
all steps voltage dependent 


c, 


o t 


AG C = -RTln^ 0/ = Q^FV, 


Noninactivating; MWC scheme; 
allosteric transition voltage 


c, 


o f 


AG C = -RT In K 0 c ^ O/ = Q n „V„ - k B T\nP™ + k B T\n(l - P 0 mi ") 


independent 








Inactivating; all steps voltage 
dependent or intermediate 


Q 


1, 


AG C = -RThyK 0 ^ = Q^FV m 


step voltage independent 








Inactivating; voltage- 
independent inactivation 


c, 


1, 


AG C = -RT In = Q m ^FV„, - RT In 7 max 



The table lists the relationship between the median measure of free energy (derived from a QV plot) and the chemical free energy difference between the 
unique initial and final states of the system AG C for six broad classes of gating modes of the ion channel. 
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Thus, Eq. 28 tells us that even in this scenario, QmaxV m is 
the change in the chemical free energy of the ensemble 
but not the chemical free energy difference between a 
unique final and initial states of the system, AG C . The 
latter will be equal to -k B Tln(K° K° ). The rela- 
tion between the two measures of free energy can be 
expressed as 



k R Tln 



1 + <r 



(29) 



The maximum fraction of inactivating channels (I max ) 
in this system, at depolarizing voltages, is / (1 + K^). 
In most inactivating ion channels, inactivation is cou- 
pled to activation, which would mean that X" » K° 0 . 
Thus, Eq. 29 can be reexpressed as 



AGc 



) V 

win ax m 



k B T In 7 max . 



(30) 



Eq. 30 is very similar to the relation that we derived for 
the case when the terminal voltage-dependent opening 
transition of a noninactivating ion channel is voltage in- 
sensitive (Eq. 23; Fig. 4 A, scheme II). The difference 
between AG C and QmaxVm is small when I max is close to 
unity and continues to increase with decreasing values 
of Imax- These points are illustrated through simulations 
performed using scheme V (Fig. 6 A) . Qma X FV m is plot- 
ted against AG C for different values of K° nact in two situ- 
ations, when inactivation is voltage dependent and 
when it is voltage independent (Fig. 6 C) . In the former 
case, they are equal, whereas in the latter case, they are 
equal only at relatively large values of K ima , where I max 
is ^T (Fig. 6 D). Thus, our overall proposition that 
QmaxFV m is the chemical free energy of the ensemble 
remains valid even in inactivating channels. 

DISCUSSION 

In summary, we have derived a general expression for 
the change in the chemical free energy change associ- 
ated with the activation of a voltage-dependent ion 
channel, AG C : 

AG C = Q mm FV m , 

where Qmax is the maximum amount of gating charge 
displacement, F is the Faraday constant, and V m is the 
median voltage of activation estimated from a QV curve. 
V m is mathematically defined by Eq. 11. It is the voltage 
at which the area bounded by the QV curve, V m axis, 
and Q = 0 axis becomes equal to that bounded by the 
QV curve, V m axis, and the Q = Qmax axis, as depicted 
pictorially in Fig. 1 B. We have considered various classes 
of models with cooperativity, voltage-independent tran- 
sitions, multiple open states, latent charge movement, 
etc. and discussed the measures of free energy change 



in each case. We demonstrate that the V m measure of 
free energy change is independent of the nature of the 
transition pathway and of the symmetry relations within 
the protein (i.e., whether parts of them are identical or 
not) . In some instances where initial or terminal steps 
of the transition pathway are voltage independent, cer- 
tain correction factors need to be added to the V m mea- 
sure to calculate the free energy difference between the 
first and last state of the channel. General relations be- 
tween the median estimates of free energy changes and 
the chemical free energy differences between the first 
and last states of the system in six broad classes of ion 
channel gating models are listed in Table 3. 

Although these relations are fairly general, their ap- 
plication to understand the energetic relationships in 
ion channels should be performed with careful delib- 
eration of the assumptions on the basis of which the 
relations are derived and the behavior of the specific 
system of interest. The median estimate of free energy 
is no substitute for full quantitative models of channel 
gating. However, it might serve as a useful check/ con- 
straint of the parameters for detail models of ion channel 
gating. Also, this measured free energy change includes 
free energy contributions from the intrinsic stabilities 
of the different structural units as well as interactions 
between them. Parsing out these individual compo- 
nents via measurement of QV curves alone is not straight- 
forward and might require the application of site-specific 
measurements along with more detailed analyses such as 
the recendy proposed x-value analysis (Chowdhury and 
Chanda, 2010). 

The median voltage of charge movement (V m ) is a 
more accurate and physically meaningful free energy cor- 
relate than Vi/2- The two parameters will be identical 
when the QV curve is symmetric, but the difference be- 
tween them can be easily appreciated in cases in which 
the QV curve is nonsymmetric and/ or shows multiple 
phases (Perozo et al., 1994; Ledwell and Aldrich, 1999; 
Tao et al., 2010; Lacroix and Bezanilla, 201 1 ) . We estimate 
that the chemical free energy change for the Shaker po- 
tassium channels is approximately —14 kcal/mol, which 
suggests that most electrical energy is being efficiently 
used to drive conformational change in the ion channel. 

Our derivations establish QV as the direct free energy 
correlate of the free energy change in a voltage-depen- 
dent system. In many instances, Po-V curves have been 
used to obtain empirical estimates of free energy of activa- 
tion of voltage-dependent ion channels (Monks et al., 
1999; Li-Smerin et al., 2000a; Hackos et al., 2002; Yrfrach 
and MacKinnon, 2002; Sadovsky and Yrfrach, 2007). Such 
energy measures are substantially different from those 
proposed in this study. Operationally, a Po-V-based esti- 
mate of the free energy change can be interpreted as the 
free energy difference between the ensemble of closed 
states of the system at the threshold voltage of pore 
opening and the final ensemble of open states. This also 
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implies that a Po-V curve does not fully capture the activa- 
tion energetics of a voltage-dependent ion channel that 
activates in multiple steps. 

On a final note, it must be mentioned that, in a case 
in which the Q mx of a channel is not measured, the 
energetic effect of a mutation on the activation ener- 
getics can be enumerated via the V m of its QV curve 
alone, under the assumption that Qmax is unaltered by 
the mutation. If the QV curve of the mutant channel 
is symmetric, V1/2 might be used instead of V m (both will 
be identical). This implies that considering the 
Boltzmann slope in the energy terms introduces a theo- 
retical anomaly in the expression for the change in free 
energy in cases beyond the realm of a two-state approxi- 
mation. Therefore when Qj„ ax is not measured, an em- 
pirical estimate of the energetic effect of a mutation on 
the voltage-gated ion channel can be assessed simply via 
comparison of the V m value of its QV curve with that of a 
control (Soler-Llavina et al., 2006; Gagnon and Bezanilla, 
2010; Muroi et al., 2010; Xu et al., 2010). The obvious 
caveats to this approach are cases in which mutations 
affect Qjnax- In such scenarios, to quantify the magnitude 
of energetic perturbation, apart from the QV curve, 
Qni ax for each charge-altering mutation needs to be 
measured separately. 
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